arXiv:1507.08136v3 [quant-ph] 27 May 2016 


Multiatom Quantum Coherences in Micromasers as Fuel for Thermal and Nonthermal 

Machines 

Ceren B. Dag,^ Wolfgang Niedenzu,^ Ozgiir E. Miistecaplioglu,^’ * and Gershon Kurizki^ 

^Physics Department, University of Miehigan, 450 Chureh St., Ann Arbor, MI fSlOO-lOfO, USA 
^Department of Chemieal Physies, Weizmann Institute of Seienee, Rehovot 7610001, Israel 
^Department of Physies, Koe University, 34450 Sariyer, Istanbul, Turkey 

In this paper we address the question: To what extent is the quantum state preparation of 
multiatom clusters (before they are injected into the microwave cavity) instrumental for determining 
not only the kind of machine we may operate but also the quantitative bounds of its performance? 

Figuratively speaking, if the multiatom cluster is the “crude oil”, the question is: Which preparation 
of the cluster is the refining process that can deliver a “gasoline” with a “specific octane”? We classify 
coherences or quantum correlations among the atoms according to their ability to serve as (i) fuel 
for nonthermal machines corresponding to atomic states whose coherences displace or squeeze the 
cavity field, as well as cause its heating; and (ii) fuel which is purely “combustible”, i.e., corresponds 
to atomic states that only allow for heat and entropy exchange with the field and can energize a 
proper heat engine. We identify highly promising multiatom states for each kind of fuel and propose 
viable experimental schemes for their implementation. 


I. INTRODUCTION 

The maser (microwave amplification by stimulated 
emission of radiation) was conceived based on thermo¬ 
dynamic considerations [1, 2]. In its micromaser imple¬ 
mentation, coherent radiation is generated by inverted 
two-level (Rydberg) atoms that are randomly injected 
into a microwave cavity one by one [3-5]. For years, the 
focus of micromaser studies had been on its quantum- 
electrodynamics features [4-16], including its extensions 
to the cooperative regime of multiatom clusters that are 
simultaneously present in the cavity [II, 12], until Scully 
et al. [17, 18] revived the interest in the thermodynam¬ 
ics of such devices. They treated the atomic beam as a 
thermodynamic resource, since randomly injected atoms, 
which are discarded (traced out) after they exit the cav¬ 
ity, constitute an effective reservoir (bath) for the cavity 
field mode (in the Markovian approximation). 

The surprising finding of Scully et al. [17, 18] was that 
a beam of three-level atoms with coherence between two 
of its levels may be viewed as a nonthermal, quantum- 
coherent (“phaseonium”) bath that, given an appropriate 
phase ip of the interlevel coherence, can thermalize the 
cavity field to a temperature > T, where T is the 
atoms’ temperature without coherence. The dramatic 
consequence of the higher temperature attainable by the 
cavity field owing to the phaseonium coherence is a trans¬ 
gression of the nominal Carnot efficiency bound in a heat 
engine, with the cavity field in the role of a working fluid 
(WF): If the WF undergoes a cycle where it is coupled to 
the phaseonium bath in one stroke and to a cold bath at 
temperature in another, then the efficiency bound of 
the engine satisfies r] < I—Tc/T^, instead of 77 < 1—Tc/T. 
This landmark proposal has triggered a variety of propos¬ 
als for engine schemes based on nonthermal baths that 
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are capable of “super-Carnot” operation [19-28], among 
them engines fuelled by a squeezed nonthermal bath [26]. 

Some of us have recently asserted [29] that machines 
fuelled by nonthermal baths may be divided into two cat¬ 
egories according to their operation paradigm: 

• Machines of the first kind are those fuelled by a 
nonthermal bath, such as a squeezed-thermal or 
coherently-displaced thermal bath, that render the 
WF steady-state nonpassive [30-36]. Such baths 
change the machine into a thermo-mechanical en¬ 
gine that, unlike a heat engine, is fuelled by both 
mechanical work and heat imparted by the bath to 
the WF. The Carnot bound may be transgressed 
in such machines at the expense of work supplied 
by the bath. However, their efficiency bound can¬ 
not be properly compared with the Carnot bound, 
since the latter is a restriction imposed by the sec¬ 
ond law on heat [37] but not on work imparted by 
the bath. 

• Machines of the second kind are those where the 
WF is thermalized by the nonthermal bath, as is 
the case of an engine fuelled by a phaseonium bath. 
Such a machine is a proper heat engine but the abil¬ 
ity of the phaseonium bath to thermalize the WF 
to a temperature > T elevates its Carnot bound 
above that associated with an incoherent bath at 
temperature T. 

Intriguingly, in micromaser setups, a beam of mul¬ 
tiatom clusters has been shown to thermalize the cavity- 
field WF in some cases [19, 20], but also coherently dis¬ 
place [21, 27] or squeeze it [38]. This implies that the 
WF may receive both work (and thus become nonpas¬ 
sive) and heat from the bath. The cavity field may thus 
be the key ingredient in machines of the two kinds sur¬ 
veyed above. However, the criteria whereby atoms in 
a micromaser can fuel machines of either the first or the 
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second kind are generally unknown, notwithstanding sev¬ 
eral recent results obtained along this line [19, 20]. 

Here we pose the question: To what extent is the quan¬ 
tum state preparation of multiatom clusters (before they 
are injected into the cavity) instrumental for determin¬ 
ing not only the kind of machine we may operate but also 
the quantitative bounds of its performance? Figuratively 
speaking, if the multiatom cluster is the “crude oil”, the 
question is: Which preparation of the cluster is the refin¬ 
ing process that can deliver a “gasoline” with a “specific 
octane”? 

To answer this question, we first derive (in Sec. II) 
a master equation that governs the cavity field under 
the standard assumption of a short interaction of each 
atom with the cavity field, ^ 1, where g is the atom- 
field coupling strength and r is the interaction time. In 
this regime, the steady-state density matrix of the cavity 
field may only be a Gaussian state: thermal, coherently- 
displaced or squeezed [39, 40]. This restriction is shown 
to imply (Sec. Ill) that two- or three-atom clusters suf¬ 
fice for the preparation of all Gaussian states of the 
field, thus making larger clusters qualitatively redundant. 
Secs. IV A and IV B are devoted, respectively, to the clas¬ 
sification of coherences or quantum correlations among 
the atoms that may be associated with (i) fuel for ma¬ 
chines of the first kind that correspond to states whose 
coherences displace or squeeze the cavity field, as well as 
cause its heating; and (ii) fuel for machines of the sec¬ 
ond kind which is purely “combustible”, i.e., corresponds 
to atomic states that only allow for heat and entropy 
exchange with the field. In both Sec. IV A and IV B 
we identify highly promising multiatom states for each 
kind of fuel and infer the best parameters relevant to 
the machine operation. In Sec. V we discuss the results 
and propose viable experimental protocols for their im¬ 
plementation. 


II. MODEL AND EFFECTIVE MASTER 
EQUATION 

We consider a micromaser-type setup wherein the 
cavity-field mode is the working fluid (WF) that is ener¬ 
gized (fuelled) by a beam of two-level atoms which are 
injected into the cavity at random, Poisson-distributed, 
times [4, 15]. By contrast to the standard micromaser 
scenario [5, 41], the atomic beam is here assumed to 
be composed of V-atom clusters that are prepared in a 
controlled quantum-correlated (entangled) state prior to 
their injection into the cavity, where they interact with 
the cavity field simultaneously (see Fig. I). Yet, the re¬ 
placement of single atoms by V-atom clusters does not 
change the basic premise of micromaser theory whereby 
their random injections allow to treat the atom clusters 
as an ergodic “bath” that is continuously coupled to the 
cavity mode, so that the latter is governed by a mas¬ 
ter equation (upon tracing out this “bath”) [5, 15, 41]. 
Nevertheless, we will show that quantum coherence or 
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Figure 1. A schematic of the two- and three-atom micro¬ 
maser model, where clusters of two-level atoms are injected 
into a single-mode cavity repeatedly in a Poissonian random 
sequence. The transition time of the atoms through the cav¬ 
ity is much shorter than the cavity lifetime, atomic relaxation 
and dephasing times or the mean free-time between the inter¬ 
actions, so that there can be at most one cluster present in 
the cavity at a time. The cavity-mode steady-state crucially 
depends on the state of the cluster, as shown here. 


interatomic quantum correlations (entanglement) in the 
cluster may crucially influence the dynamics of the cav¬ 
ity field. Experimentally feasible schemes for the present 
scenario will be discussed in Sec. V. 

In keeping with the standard assumptions of micro¬ 
maser theory [16, 41], we take the transition time of the 
atoms through the cavity to be short enough to neglect 
atomic relaxation and dephasing as well as cavity loss, 
and to assume that there can be at most one cluster 
present in the cavity at a given time. Under these stan¬ 
dard assumptions, we may derive a master equation for 
the dynamics of the cavity field. 

The interaction of the atomic cluster with the cavity 
is described by the Tavis-Cummings model [42] 


+ ( 1 ) 

where the atom, cavity and interaction Hamiltonians are 
respectively given by 



II 

(2a) 
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He 

= huJeCi^a; 

(2b) 
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Here a, are the annihilation and creation operators for 
the cavity field and (j|, aj' are the z, raising and low¬ 
ering Pauli operators for the jth atom with j = I,..., Y. 
The atomic transition frequency cja is resonant with the 
cavity frequency ujc- The interaction between the atoms 
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and the cavity is assumed to be spatially homogeneous 
with strength g. 

Under the foregoing assumptions, the combined sys¬ 
tem of the atomic cluster and the cavity field evolves 
unitarily during the short interaction time r. The uni¬ 
tary propagator U{r) = exp(—in the interaction 
picture can be analytically computed to second order in 
gr (see Appendices A-C). Denoting the injection time 
of the jth cluster into the cavity by tj, the evolution of 
the reduced density operator of the field mode, which is 
obtained upon tracing out the atoms, reads [43, 44] 

Pitj+T) = T^a[UiT)pa® p{tj)U'’{T)] = S{T)p{tj). (3) 

Here pa is the initial density operator of the atomic clus¬ 
ter and S{t) is a superoperator that propagates the cav¬ 
ity state p{tj) to p{tj + r). The atomic clusters arrive 
randomly at a rate p and pass through the cavity within 
a time interval of with a probability of pSt. The 

field changes according to S{t) when a cluster is present 
and otherwise does not change at all, so that the overall 
change of the field state is 

-f St) = pStS{r)p{t) + (1 - pSt)p{t). (4) 

For St ^ 0 we obtain the master equation 

p{t) =p[S{t) -l\p{t) (5) 

which describes the Markovian dynamics of the single¬ 
mode cavity [5, 20-22]. Here we do not include the usual 
cavity decay term, in order to clearly identify the role of 
coherences for the field evolution, particularly whether it 
thermalizes or not, but it is straightforward to do so. 

The master equation (5) can be rewritten as 


p{t) = p 


N 


N 


UniiT)pit)[UnjiTp - pit) 




n=l 


, (6) 


where denote the matrix elements of pa- In the stan¬ 
dard basis of energy-state products the diagonal elements 
an are the populations and the off-diagonal elements aij 
with i j are coherences or correlations, respectively. 
This standard basis is {|e) , |p)} for one atom and shown 
in Fig. 2 for two ({|ee), \eg ), \ge ), \gg)}) and three atoms, 
respectively. 

Using the explicit forms of the respective propa¬ 
gators U{t) for one-, two- and three-atom clusters 
[Eqs. (A2), (B9) and (C6) in the Appendix], the master 
equation (6) can be expressed in the illuminating form 

p^ -i [Hes, p] + Kp + hp. (7) 

Here the first term corresponds to the effect of a coherent 
drive applied to the cavity, which is described by the 
effective Hamiltonian 


(8) 




Figure 2. Energy levels of clusters of two- and three two level 
atoms. The numbers next to the levels correspond to the 
indices used in the text to denote their corresponding position 
in the natural basis. 


The Lindbladian in Eq. (7) describes a squeezing pro¬ 


cess and is given by 

w = (9) 

where p = p{gr)‘^ is an effective coupling rate. The 
squeezing excitation and de-excitation Lindbladians are 

= 2a^pa^ — p — pa^a^ (10a) 

Lf = 2apa — aap — paa^ (10b) 

respectively [39]. The Lindbladian L is given by 

hp = p (^Le + , (11) 

where 

= 2apa^ — a^ap — pa^a (12a) 

Lg = 2o) pa — aa^ p — pao) (12b) 


are the Lindbladians for incoherent de-excitation and ex¬ 
citation, respectively. The coefficients for different clus¬ 
ter sizes are shown in Table 1. 

The master equation (7) allows for the generation of ar¬ 
bitrary (Gaussian) field states, i.e., thermal-, displaced- 
and squeezed states. Higher-order (i.e., non-Gaussian) 
processes cannot be induced by a second-order mas¬ 
ter equation. The case A" = 2 is the minimal clus¬ 
ter generating these processes: Adding another particle 
{N = 3) does not make a qualitative difference compared 
to Eq. (7). 


III. CLASSIFICATION OF COHERENCES AS 
DIFFERENT TYPES OF FUEL 

The key observation we infer from Eq. (7) is that co¬ 
herences or correlations in the multiatom cluster may be 
classified according to the disjoint blocks in the density 
matrix p^ that are associated with qualitatively different 
terms in the master equation, each term giving rise to a 
different kind of field dynamics. Eiguratively, the differ¬ 
ent coherences are different types of fuel for the cavity- 
field WE (Eig. 3): 


He^i = pgr (Aa'^ + \*a) . 



















4 



1 atom 

2 atoms 

3 atoms 

Ve 

ail 

2aii + 022 + ^33 + a23 + ^32 

3aii + 2 De + Dw + Ce + Cw 

^9 

022 

2a44 + 022 + O33 + a23 + <^32 

3o88 + 2Dw + De + Ce + Cw 

X 

ai2 

ai2 + ai3 + a24 + 0-34 

025 + 033 + a46 + a47 + a26 + 037 
+ai2 + ai3 + ai4 + + oqs + 073 


0 

Ol4 

028 + <^38 + a48 + Oi3 + ttie + O17 


Table L Coefficients of the master equation (7) for different cluster sizes. For later convenience we have defined for three- 
atom clusters the abbreviations De = CL22 + ^33 + 0^44, Dw = ^55 + ^66 + ^rr, Ce — CL23 + ^24 + <^32 + <^34 + ^42 + 0^43 and 
Cw — OLbQ + 0^65 + ^57 + ^75 + CiQ7 + CL7Q- 


• The blocks adjacent to the main diagonal of the 
Pa matrix in the standard basis of Fig. 2 contain 
coherences that can only induce absorption- and 
emission processes in the field (WF), as they are as¬ 
sociated to Lp in the master equation (7). We shall 
refer to these elements as heat-exchange coherences. 
They have a caloric (“flammable”) value, i.e., they 
may contribute to the thermalization of the cav¬ 
ity field. Heat-exchange coherences do not arise in 
the single-atom case, as they correlate states of the 
same energy, e.g., \eg){ge\ in two-atom clusters and 
|eep)(epe| in three-atom clusters. 


I Ineffective coherences 

1 Populations 

S 

2 Heat-exchange coherences 
Displacement coherences 


three atoms 



• Displacement coherences associated with the 
—i[F7eff, p] term in the master equation (7) arise for 
all cluster sizes as they correlate states differing by 
one excitation, i.e., |e)(p| and |p)(e| in single atoms, 
\eg){ge\ and its Hermitian conjugate in two-atom 
clusters and, say, \eeg){geg\ in three-atom clusters. 

• Squeezing coherences correspond to an exchange of 
two excitations and may exist in two-atom clusters 
in the form of \ee){gg\ and its Hermitian conju¬ 
gate, or in three-atom clusters in, say, the form 
\eeg){ggg\. 

Those matrix elements of pa that do not contribute to 
the field evolution shall be called ineffective coherences. 

Larger cluster sizes will not change the qualitative fea¬ 
tures of the master equation (7), which holds to second 
order in gr and thus may only induce the same second- 
order (Gaussian) processes as listed above. The different 
types of coherences and their relation to the number of 
excitations is illustrated in the tree diagram of Fig. 4 
(which should not be confused with the same term in 
graph theory). 


IV. CORRELATED ATOMIC CLUSTERS AS 
FUEL FOR MACHINES OF THE FIRST AND 
SECOND KIND 


I Squeezing coherences 


Figure 3. Density matrix of the atomic cluster for one (left), 
two (middle) and three (right) atoms, respectively, with color- 
and pattern-filled squares representing the different roles of 
the coherences with respect to the cavity-field evolution de¬ 
scribed by the master equation (7). Red plain dark squares 
are populations and light blue squares are ineffective coher¬ 
ences. Yellow diagonal striped squares are zeroth order co¬ 
herences that can contribute to thermalization. Dark blue 
vertical striped squares are first order coherences that can 
contribute to the coherent displacement of the cavity field. 
Green horizontal striped squares are second order coherences 
contributing to the squeezing of the cavity field. 


that not only heat but also work has been trans¬ 
ferred from the bath to the cavity mode. Con¬ 
sequently, a machine fuelled by such a bath is a 
machine of the first kind that operates thermo- 
mechanically. 

• If the atomic state only contains heat-exchange co¬ 
herences, the mode is thermalized by the bath and 
only heat is exchanged. Such a setup is thus a vi¬ 
able implementation of a heat engine powered by 
a nonthermal bath, which has been dubbed a ma¬ 
chine of the second kind. 


The beam of atomic clusters interacting with a cavity 
mode can realize one of two operation paradigms [29]: 

• If displacing or squeezing coherences are present 
in the bath, the cavity state becomes nonpassive 
(displaced or squeezed, respectively), which implies 


A. Conditions for fuelling machines of the first kind 

Under what conditions does the master equation (7) 
possess a nonpassive (nonthermal) steady state of the 
cavity mode that is required for a machine of the first 
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Figure 4. Trees of coherence for one- (left), two- (middle) 
and three- (right) atom clusters, respectively. The circles are 
the basis states, where same colour indicates the same num¬ 
ber of excitations. Blue solid lines indicate heat-exchange co¬ 
herences that may contribute to thermalization of the cavity 
field. These coherences between states with the same number 
of excitations only appear in the multipartite case. Displace¬ 
ment coherences (red dotted lines) between states differing by 
one excitation arise for all particle numbers. Squeezing coher¬ 
ences (green dashed lines) are between states differing by two 
excitations. 

kind (a thermo-mechanical machine)? What is the na¬ 
ture of such a state? 

To obtain an insight into these questions, we write the 
(Ehrenfest) equations of motions of the field mean value, 
variance and mean intensity (photon number) in terms 
of the Pa matrix elements that are grouped into the co¬ 
efficients r^, Te, A and ^ as detailed in Table I, 

{d{t)) = -^{rg- Te) {a{t)) - ipgrX, (13a) 

(d^(l)) = -p{rg - Te) {a^{t)) - 2ipgT\ {a{t)) - 2/i^, 

(13b) 

(n(i)) = -pL{rg - re) {n{t)) + pre 

- ipgT{\ {a\t)) - A* {a{t))), (13c) 


whose steady-state solutions read 



We see that a nonzero A increases the thermal mean pho¬ 
ton number by coherently displacing the cavity field to a 



Figure 5. A doubly-excited state (top) gives rise to a squeezed 
state of the cavity field. By contrast, a triply-excited state 
(bottom) thermalizes the cavity to an ultrahigh temperature. 

nonzero expectation value [Eqs.(14a) and (14c)]. Accord¬ 
ingly, the cavity field attains thermal-coherent character. 
By contrast, a nonzero ^ introduces quadrature squeez¬ 
ing to the cavity field and hence the cavity field acquires 
a squeezed-thermal character. 

As an example, consider the two-atom state 

= cos'?^ Ipp) + sin'i9 |ee). (15) 

Eor sin^ ^ ^ (for simplicity we restrict the angle to 

0 < '^9 < 7r/4) this state gives rise, according to Eq. (7), 
to the master equation (cf. Eig. 5) 

p = pi sin(2'i9)Lsp + p [cos^ -\- sin^ ^ Lg] , (16) 

which may be cast into the standard form that yields 
thermal-squeezed solutions [39, 45] 

p = nMl^sP + ^ [(^ + 1)L^ + Nhe ], (17) 

with the coefficients 

N = n(cosh^ r + sinh^ r) + sinh^ r (18a) 
M = coshr sinhr(2n + 1). (18b) 

Here r denotes the squeezing parameter and n the am¬ 
bient photon number. Upon comparing Eqs. (16), (17) 
and (18) we find 


K = pcos(2'i9) 

(19a) 

n = 0 

(19b) 

r = atanh (tan P ). 

(19c) 


In conventional experimental squeezing schemes, squeez¬ 
ing parameters range from r ^ 0.4 [46, 47] up to r ~ 
1.46 [48]. Remarkably, the squeezing parameter (19c) 
may greatly surpass existing values if we choose P ^ tt/4 
(see Eig. 6). 

Let us note that the existence of a steady-state solu¬ 
tion (whether coherent or not) requires Vg > rg, implying 
sin^ P < ^. We will discuss this point in more detail in 
what follows. 


B. Conditions for fuelling machines of the second 
kind 

It is thus clear that thermal equilibrium of the cav¬ 
ity field requires A = ^ = 0 so that all coherent pro¬ 
cesses vanish. This can be achieved by either setting 
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the displacement- and squeezing coherences to zero in 
Fig. 3 or making their respective contributions to the co¬ 
efficients A and ^ (see Table I) to cancel each other. Un¬ 
der these conditions that are required for thermalization 
of the cavity field, the master equation reads 

( 20 ) 

where /ir^/2 is the rate of emission of quanta into the 
bath and /ire/2 is the absorption rate of quanta from the 
bath. By virtue of the Kubo-Martin-Schwinger (KMS) 
detailed-balance condition [45], it is then possible to at¬ 
tribute a temperature T to the effective bath through 


Figure 7. Steady state mean number of photons in a cavity 
pumped randomly with three atom clusters in W class states, 
parameterized with angular variables 0 and 'ijj, when (5 = 0 
and (/) = 0. The symmetric W state yields the largest mean 
photon number in equilibrium and hence can be imagined as 
the "hottest" effective three atom reservoir among the W class 
states. 


this steady-state solution only exists if 0 < T < oo, 
i.e., if Tg > Te- The case Vg = re, formally resulting in an 
infinite bath temperature, can be identified as the maser 
threshold [49] (see also Appendix D). 


This temperature T is only positive and finite if rg > re. 

Let us pretend that this effective bath is composed of 
fictitious oscillators. Then Eq. (21) yields 


where 


_ n 

Tg n + 1 ’ 



( 22 ) 


(23) 


is the effective thermal excitation of the bath (mean num¬ 
ber of fictitious quanta) at the cavity frequency ujc> The 
denominator in Eq. (22) then corresponds to the sum of 
stimulated and spontaneous emission into the effective 
heat bath. 

This effective description of the effective bath is vin¬ 
dicated by the steady-state solution of the master equa¬ 
tion (20) which is the thermal (Gibbs) state [45] 


= (24) 

of the cavity mode, where Z denotes the partition func¬ 
tion. Due to the unbounded character of the Hamiltonian 


We have thus arrived at an important conclusion: A 
nonthermal beam of atoms interacting with a single cav¬ 
ity mode may act as an effective heat bath for the lat¬ 
ter, thereby thermalizing it to a finite temperature T, al¬ 
though the quantum state of the atoms may be distinctly 
nonthermal^ i.e., the atomic-cluster state is not associated 
with the notion of temperature. Nevertheless, it will drive 
the cavity field mode into a Gibbs state with a finite and 
positive temperature provided the cavity mode is below 
the maser threshold. This conclusion is consistent with 
the well-known fact that the regime below the micro¬ 
maser threshold is thermal radiation with a thermody¬ 
namic equilibrium temperature [49, 50]. Here, however, 
this temperature T of the cavity field depends explicitly 
on the coherences and correlations of the atoms. 


In Table H we present the explicit dependence of 
the temperature T, the steady-state photon number 
(n)ss, and the micromaser threshold on the multiatomic 
density-matrix parameters from Table I. 


In what follows, we focus on the relation between the 
correlations in distinctly entangled states of the cluster 
on the temperature and threshold conditions. 
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1 atom 

2 atoms 

3 atoms 

huc j \ 


In ( ^“ 44 -|-a 22 -|-a 33 -|-a 23 -|-a 32 5 1 
\ 2011+022+033-1-023-1-032 yj 

Tin f 3o88+2+>w++>£;+C'\ 1 

[ \ 3oii+2D£;+Di 4/+C' y 


“22 

2044 + 022 + 033 + 023+032 

3 o 11+ 2 +)£;++>14/+C 

Wss - 

022—0-11 

2(044 — 011) 

3(o88 —oii) + +)vV^— 

valid for {rg > re) 

Ci22 > CLii 

a44 > ail 

3 a 88 + Dw > 3 aii + De 


Table II. Steady-state properties following from the thermal master equation ( 20 ) for different cluster sizes. Here we have 
defined C — Ce + Cw- 


1 . Cavity thermalization via singly-excited entangled 
three-atom states 


We may parameterize the singly-exited entangled 
states of three atoms via 

l^)gen = COS 0COS ^ \egg) + sin6» cos \geg) 

+ sin \gge). (25) 

The mean photon number in terms of the angles then 
reads (cf. Table II and note that Ce and De vanish) 


2 


(n)gg = 1 + sin 20 cos^ ^ cos 0 + cos 0 sin 'ip cos 6 

+ sin^ sinT/^ cos(0 — J). (26) 


The maximum photon number corresponds to ^ = 0 and 
0 = 0. The variation of the mean number of photons with 
respect to the remaining parameters 0 and 0 is shown in 
Fig. 7. It is seen that the maximum number corresponds 
to the symmetric W state 

W) = {\gge) + \geg) + \egg)) /\/3, (27) 

which is known for its robust entanglement [51]. This 
state yields = 3, which means that coherences (cor¬ 
relations) in the W state increase the photon number in 
thermal equilibrium from the value {n)^ss = 1, which 
would be the case for a phase-aver aged W state. 

This amplification of photon population in the cav¬ 
ity is here due to Dicke superradiance [52]: The quan¬ 
tum interference in the W state enhances the processes 
described in the master equation, which equilibrate the 
cavity field to a canonical thermal state. Although each 
cluster is in a pure state, the entropy of the cavity in¬ 
creases via the partial-trace operation after each interac¬ 
tion, which removes the information about the atomic 
state. The crucial contribution of the heat-exchange 
coherences present in the W state (27) can be traced 
to the effective temperature T ^ 3.47/ico’//cb that these 
coherences induce as to compared to the temperature 
obtained for its phase-averaged (classically-correlated) 
counterpart (i.e., without any heat-exchange coherences). 
To ~ lAAhuj/kB < T. This temperature Tq is solely de¬ 
termined by the populations of the computational-basis 



TT Stt Stt Stt 

2 4^42 
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Figure 8 . Mean photon number ( 29 ) in a cavity pumped by 
atom clusters in the generalized GHZ state ( 28 ). 


states and may hence be thought of a “classical” effect. 
By contrast, the augmented temperature T > Tq stems 
from the heat-exchange coherences (that here lead to con¬ 
structive quantum interference). The deviation from Tq 
is thus of quantum-mechanical origin. 

Therefore, we conclude that the symmetric W state 
provides the highest equilibrium temperature to the cav¬ 
ity field among the entangled singly-excited states. This 
comes about since in the symmetric W state all contri¬ 
butions i^int |^)gen Coherently, allow¬ 

ing for cooperatively enhanced interaction in this three- 
particle Dicke state [53]. 


2 . GHZ states: Towards infinite effective temperature 


Equation (13c) only possesses a steady-state solution 
[Eq. (14c)] if ass > (cf. Table II). This condition can 
be fulfilled for a generalized GHZ state parameterized by 


|GHZ)g^„ = cos ^ |eee) + sin ^ \ggg ), (28) 

for which the mean photon number in thermal equilib¬ 
rium becomes 


We have plotted the steady-state photon number as a 
function of 'P in Eig. 8. The figure shows that as 'P ^ 7r/2 





























one reaches the micromaser threshold, where the mean 
photon number diverges, signifying an infinite effective 
temperature of the bath. 

Let us now consider the case'd = 7r/2 that yields the 
GHZ state 

|GHZ) = (| 555 ) + |eee))/V2. (30) 

The only nonzero coherences are then aig and agi, which 
according to Fig. 3 are ineffective. Indeed, the only 
nonvanishing parameters of the master equation (7) are 
= 3/2. Such parameters in the thermal Lindblad 
equation 

= ^(Le+Ld) (31) 

correspond, according to the KMS detailed-balance con¬ 
dition, to an infinite temperature of the bath. A beam 
of GHZ states thus cannot thermalize a cavity mode. In¬ 
deed, the mean photon number in the cavity grows in an 
unbounded fashion, 

{n) = f, (32) 

according to Eq. (13c). 

Similarly to the symmetric W state, a nearly symmet¬ 
ric GHZ-type state (cf. Fig. 5) is the optimal choice for 
reaching high cavity temperatures. However, the mecha¬ 
nism is entirely different in the two cases. The symmetry 
in the W state allows for constructive quantum interfer¬ 
ence (superradiance), so that the enhancement is purely 
quantum-mechanical. By contrast, the nearly-symmetric 
entangled GHZ state allows us to approach the maser 
threshold. The photon number (29) only depends on the 
populations since, according to Fig. 3, all coherences in 
the state (28) are ineffective. Hence, the phase-averaged 
counterpart of Eq. (28) results in the same effective tem¬ 
perature as that of the state (28), i.e., Tq = T. 

3. States leading to ultrahigh temperatures of the eavity field 

Let us incoherently mix W states with states that we 
denote as E states that belong to the general class of 
three-atom W states with two excitations in the upper 
blue triangle in Fig. 4 (which are also Dicke states [54]), 

|E) = T {\eeg) + \ege) + \gee)). (33) 

According to Table H they contribute to Ce and De and 
correspond to Tg = 4 and = 3, namely their rate of 
absorption surpasses the emission rate, leading [by the 
KMS detailed balance condition (21)] to a negative tem¬ 
perature^ which is outside the scope of this paper. 

The chosen nearly equal mixture of W- and E-states 
has the form 

PWE= (^l+eyW){W\ + (t-eyE){E\, (34) 


where 0 < 5 <C 1 is a small positive number. This state 
is a mixture of the two coherence triangles indicated by 
solid blue lines in Fig. 4. The corresponding nonvanish¬ 
ing parameters of the master equation (7) rg = 7/2 — 5 
and rg = 7/2 + 5 (cf. Table I) imply that > rg, so 
that this mixed state corresponds to a positive and finite 
effective temperature, T « lehuj/AkB^ The correspond¬ 
ing classical-like phase-averaged counterpart of Eq. (34) 
would thermalize the cavity to Tq = 3shuj/4kB- Here, 
two enhancement factors are involved. The first factor 
is C (cf. Table H), enhancing Tq to higher temperatures 
by the quantum interferences due to the coherences in 
the W and E states. The second factor is the classical 
enhancement of Tq due to the operation near the maser 
threshold. 


V. DISCUSSION 

We have studied the thermodynamic implications of a 
generalized micromaser model wherein the cavity mode 
interacts with a beam of quantum-coherent or quantum- 
correlated multiatom clusters. Our central goal has been 
to classify the states of such clusters prior to their in¬ 
jection into the cavity according to their ability to fuel 
the cavity field as “working fluid” in a machine of either 
the first kind (thermo-mechanical engine) or the second 
kind (heat engine). To this end we have derived a Lind¬ 
blad master equation for the cavity field mode that de¬ 
scribes absorption- and emission of the field, its coherent 
displacement and squeezing caused by the atoms that 
may act, respectively, as a thermal, displaced-thermal 
or squeezed-thermal bath. These distinct Gaussian pro¬ 
cesses that the field may undergo are determined by the 
prefactors of the respective terms in the master equation 
that are, in turn, determined by disjoint blocks (coher¬ 
ences) of the multiatom density matrix. 

The main results of our analysis are as follows: 

• An important insight that we have obtained is that 
two- and three-atom clusters are capable of acting 
as fuel for both kinds of machines in a highly ef¬ 
fective fashion, so that there is no need to involve 
larger clusters. Still, a larger number of coherences 
as the cluster grows in size may further enhance the 
work output. 

• For machines of the first kind, our analysis has re¬ 
vealed a particularly promising, simple, fuel in the 
form of two-atom clusters whose state is a nearly 
equal superposition of doubly-excited and doubly- 
ground states. Such a state is expected to give 
rise to very large squeezing of the cavity field. It 
may thus present a far superior alternative to ex¬ 
isting squeezing schemes of cavity fields [46-48]. 
Such a strong squeezing may have fascinating ap¬ 
plications [48] also outside of quantum thermo¬ 
dynamics. Our interest here is that this strong 
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squeezing source may fuel a cavity field in a hy¬ 
brid thermo-mechanical machine [29] with nearly 
100% efficiency, at the expense of mechanical work 
supplied by the two-atom clusters. 

• For machines of the second kind, we have found W 
states of three-atom clusters to act as conventional 
heat-baths fuel at a positive finite temperature that 
is controllable by the W state. By contrast, three- 
atom GHZ- and E-states have been found to cor¬ 
respond to effective baths at infinite or negative 
temperatures, respectively, that do not allow for a 
thermal steady-state solution for the cavity field. 
On the other hand, nearly-equal mixtures of W 
and E states have been identified as fuel capable 
of thermalizing the cavity field to an ultrahigh tem¬ 
perature. 

To conclude, our results are potentially useful for the 
design of thermal and nonthermal machines based on 
micromaser setups. The availability of all Gaussian 
processes via preparation of two- and three-atom clus¬ 
ters allows to implement heat engines (wherein the cav¬ 
ity field is thermalized) but also thermo-mechanical en¬ 
gines (wherein the cavity mode is coherently displaced or 
squeezed). 

We wish to stress the feasibility of the diverse forms 
of state preparation of multiatom clusters (prior to their 
injection into the cavity) we have employed in our anal¬ 
ysis: 

• The arsenal of quantum gate operations [51] can in 
principle prepare two or three trapped atoms in an 
entangled state on demand, but such preparation 
may require single-atom addressability. 

• Alternatively, W states can be generated via quan¬ 


tum feedback control [55] or at fusion-based light- 
matter interfaces [56]. Multipartite entangled 
states may also be generated via photon-mediated 
interactions, as recently discussed in [57]. 

• Another alternative is an optimized probabilistic 
scheme for multiatom entangled-state preparation 
in a cavity [58]. 

• Eor two-atom entangled-state preparation we may 
resort to controlled diatomic dissociation [59], col¬ 
lisions in a cavity [60] or long-range dipole-dipole 
interactions [61]. 

On the fundamental side, our results provide clues to 
the thermalization or nonthermalization of a system (here 
the cavity field) via its contact with quantum-correlated 
multipartite clusters that act as nonthermal baths. Such 
processes reflect the subtle rapport between quantum 
correlations in the bath and thermalization [20]. 
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Appendix A: Time-evolution operator for a one-atom micromaser 

Eor one atom the time-evolution operator 

U{t) = exp(-iiJintr) 

to second order in gr readily evaluates to 

0 a\ _ A ~ + 1) —igra 




2 0 


-igra 


^ 1 - 


(Ai) 


(A2) 


Appendix B: Time-evolution operator for a two-atom micromaser 


The time-evolution operator 


U{t) = exp(-iiJint'r) = exp(-i5TP) 


(Bl) 
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of the joint cavity-atoms system can be computed to second order in gr using the collective angular-momentum 
operators S± = Y^j=i o'f such that 

P = a5'++a'^5'_. (B2) 

The latter can be decomposed into irreducible subspaces by changing from the computational basis (spanned by 
products of single-atom |e) and \g) states, cf. Fig. 2) to the basis of Dicke states [53] by means of the transformation 
matrix 


T = 


/ 1 0 0 0 \ 

Vo 0 1 0 / 


The addition of two spin- 1/2 gives rise to a triplet and a singlet, 505 = 1 © 0. As a consequence, 

/ 0 V2a 0 0\ 

T^PT = 


\/2at 0 V2a 0 

0 \/ 2 at 0 0 

VO 0 00 / 


Pi 


and the propagator is given by the direct sum 

where to second order in gr 


U{t) = Ui{t)(BUo{t), 

(gr? 


Ukir) « Ijk - igrPk - 

Here Ik denotes the unit matrix of the same dimensionality as Pk. Explicitly, we find 

1 


Ui{t) = 


(gr^aU + l) 
—i\/2a^gT 1 

-a^^igrf 


—i\/ 2 agT —o?{gT)‘^ 

(^r)^( 2 a’t'a + 1 ) —i^/2agT 
—i\/2a^ gr 1 — {grYa^a 


and 


t/o(r) = 1. 

Transforming back to the computational basis {|ee), \eg ), \ge) , \gg)} yields 


U{t) = 


/I — {gT)‘^{o)a + 1 ) —iagr —iagr 

1 - 5(5'r)^(2a'fo + 1) -\{gTY{2a^a + 1) 
1 ) I - ligrf {2a)a + 1) 


—la^gr 
—ia^gr 


V 




aPigrf 


—ia^gr 


—ia)gT 


-a^grf \ 
—iagr 
—iagr 

1 - {gTfo)a J 


(B3) 

(B4) 

(B5) 

(B 6 ) 

(B7) 

(B 8 ) 

(B9) 


Appendix C: Time-evolution operator for a three-atom micromaser 


For three particles one proceeds exactly like in the preceding section. The transformation matrix now reads [53] 



1 

0 

0 

0 

0 

0 

0 

0 

\ 


0 

1 

x/3 

0 

0 

0 

0 

7 

0 



0 

1 

73 

0 

0 

1 

72 

0 

76 

0 



0 

1 

73 

0 

0 

1 

72 

0 

1 

76 

0 



0 

0 

1 

©3 

0 

0 

1 

/2 

0 

1 

©6 



0 

0 

1 

73 

0 

0 

1 

72 

0 

1 

76 



0 

0 

1 

73 

0 

0 

0 

0 

[2 

V 3 



0 

0 

0 

1 

0 

0 

0 

0 

/ 


(Cl) 











11 


The addition of three spin-1/2 gives rise to a quadruplet and two doublets, ^ 


| = |0|0|. As a consequence, 


T^PT : 


/ 0 

\/3a 

0 

0 

0 

0 

0 

0 \ 

V^at 

0 

2a 

0 

0 

0 

0 

0 

0 

2at 

0 

\/3a 

0 

0 

0 

0 

0 

0 

^/3a^ 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 

0 

0 

0 

0 

0 

at 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 

V 0 

0 

0 

0 

0 

0 

at 

0/ 


3/2 


1/2 


(C2) 


1/2 


and the propagator is given by the direct sum 

U{r) = Us/ 2 {r) 0 Up^ir) 0 t/i/ 2 (r), 

where 


(C3) 


UzMt) = 


/ 1 - 1 ) 

—i\/3a^gT 

0 


—i\/3agT 

1 - 5(fl'r)^(7a^a + 4) 
—2ia)gT 
-\/3a^^(5r)2 


-V^a'^igTf 

—2iagT 


0 \ 

-\/ia^{gTf 


1 — |(gr)^(7a^a + 3) —i\/3agT 

—i^/Sa^gr 1 — ^{gT)^a^a j 


(C4) 


and 


1 — ^(^r)^(a’*'a 0 1) —iagr 

0r)^a0 
2 


—ia^gr 


1 - 


C7i/2(r) = 

Transforming back to the computational basis (cf. Fig. 2) yields the matrix elements 


(C5) 


C7ii = ^ (2 - 3(5T)2(a^a + 1)) , 

U 21 = Usi = 1/52 = Uq2 = Uss = 1/73 = F 41 

= Us5 = ^64 = b^74 = Use = %7 = —igTo^^ 

Ui2 = t/l3 = U 25 = b^26 = Use = b^37 = ^14 

= Ues = ^46 = K 47 = t/es = ^78 = —icigT 

U51 = Uei = t/71 = Us2 = Us 3 = Us 4 = -{grf 

Ui 3 = t/i 6 = Uu = U28 = U38 = U48 = -{grfa^ 

U 22 = U 33 = t /44 = 1 - ^{grfiSa^a + 2 ), 

Us 2 = U 42 = U 2 S = b^43 = Uee = Uje = b^24 = UsA 
= Use = Uje = Uer = ^67 = ~2^9^)‘^{2a^a 0 1), 
7/55 = 7766 = 7/77 = 1 - i(5T)2(3ata + 1), 

7788 = 1 - ^(5'r)^a’''a 

of the time-evolution operator. The remaining elements evaluate to zero. 

Appendix D: Maser threshold 


(C6) 


r 


It is illuminating to derive the threshold condition from 
the mean photon number, following Ref. [20]. We first 
take A = ^ = 0 and thereby eliminate the coherent- 


displacement and squeezing terms. According to Eq. (3), 
the cavity density matrix will change to 


p{tj + t) « {grf + p{tj) (Dl) 
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after the passage of the jth atom during the short interac¬ 
tion time r. The mean photon number can be calculated 
to be 

{n{tj + r)) = Tr [p{tj + T)n] = k + {gT^^. 

(D2) 

The change of the mean photon number between consec¬ 
utive injections of two atom clusters is determined by the 
increment ratio k, which is given by 

k = l-{gTf'^^. (D3) 


Assuming that the cavity is initially in the vacuum state, 
the last term in Eq. (D2) yields the mean number of pho¬ 
tons after the first-cluster passage, {n{r)) = (^r)^re/2. 
After the jth cluster passage the mean photon number 
rises to 

(n(t,)) = ^fc*-qn(r)). (D4) 

i=l 

The summation in Eq. (D4) is convergent if /c < 1, which 
is equivalent to the threshold condition Vg > Vq. As j 
oc, the summation converges to Eq. (14c). 
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